Draft version June 2, 2009 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



A THIRD GIANT PLANET ORBITING HIP 14810 
J. T. Wright 

Department of Astronomy, 226 Space Sciences Building, Cornell University, Ithaca, NY 14853 

jtwright@astro.cornell.edu 



o 
o 

(N 
(N 

6 



> 

in 
o 

O 

OS 

o 



X 



D. A. Fischer 

Department of Physics and Astronomy, San Francisco State University, San Francisco, CA 94132 

Eric B. Ford, D. Veras, J. Wang 

Department of Astronomy, University of Florida, 211 Bryant Space Science Center, P.O. Box 112055, Gainesville, FL 32611-2055 



G. W. Henry 

Center of Excellence in Information Systems, Tennessee State University, 3500 John A. Merritt Blvd.. 



Box 9501, Nashville, TN 37209 



G. W. Marcy, A. W. Howard 1 

Department of Astronomy, 601 Campbell Hall, University of California, Berkeley, CA 94720-3411 

AND 

John Asher Johnson 2 

Institute for Astronomy, University of Hawai'i, Honolulu, HI 96822 
Draft version June 2, 2009 

ABSTRACT 

We present new precision radial velocities and a three-planet Keplerian orbit fit for the V = 8.5, 
G5 V star HIP 14810. We began ob serving this star at Keck Observatory as part of the N2K Planet 
Search Project. IWright et alJ ( 2007f > announced the inner two planets to this system, and subsequent 
observations have revealed the outer planet planet and the proper orbital solution for the middle 
planet. The planets have minimum masses of 3.9, 1.3, and 0.6 Mj U p and orbital periods of 6.67, 147.7, 

and 952 d, respectively. We have numerically integrated the family of orbital solutions consistent with 
the data and find that they are stable for at least 10 6 yr. Our photometric search shows that the 
inner planet does not transit. 

Subject headings: planetary systems 
1. INTRODUCTION 



stars: individual (HIP 14810) 



The first multiple exoplanet system detected around 
a normal star 3 wa s the triple system v Andromadae 
(jButler et al.l fl999). Today, over 3 systems comprisin g 
more than one planet are known ([Wright et alJ T2009). 
including seven triple syst ems, two quadruple systems 
(fi Arae (|Pepe et al.1 120071 ) and GJ 581 (Mayor et al. 
AfcA submitted)) a nd the quintuple system 55 Cancri 
(jFischer et al.ll2008ft . 

Individual multiplanet systems offer insights into the 
dynamical evolution of planetar y systems that single- 
ton sy stems cannot. For instance. iFord. Lvstad. k, Rasiol 
(2005) showed that v Andromedae bears the scars of 
strong planet-planet scattering events preserved in their 
planets' orbital parameters. Other systems show ev- 
idence of migration and eccentricity pumping through 
mean-motion resonances (MMRs), which may be the 
sign post of convergent migration i n multi-planet systems 
(e.g. Kiev. Peitz fc Brvdenll2004D . Multiplanet systems 

1 Townes Fellow, Space Sciences Laboratory 

2 NSF Postdo ctoral Fellow 

3 Prior to this, Wolszczan & Frail ( 1992), detected the first plan- 
ets outside the solar system: three extraordinary planets orbiting 
the pulsar PSR 1257+12. 



with planet-planet interactions strong eno ugh to be de- 
tected at current RV precision (as in GJ 876. lRivera et al.l 
|2005|) allow for measurement of the inclination of the sys- 
tem, providing true planet masses. 

Comparison of multiplanet systems as an ensemble to 
singleton systems provides observational constraints to 
theories and models of the early dynamical evolution and 
migra tion history of planetary systems. IWright et al.l 
(2009) showed that the while the eccentricity distribu- 
tion of planets in multi-planet systems is similar to that 
of apparently singleton systems, their semimajor axis dis- 
tributions differ significantly. The concentration of plan- 
ets with orbital periods near 3-days seen in the single- 
planet systems is absent in multi-planet systems, as is 
the sharp jump beyond 1 AU in planet frequency. The 
fact that these features (the 3-day pileup and the 1 AU 
jump) are functions planetary multiplicity strongly sug- 
gests that planet-planet interactions play a key role in 
both migration and the origin of eccentricities. 



2. THE TENTH TRIPLE SYSTEM 

The announcement of HIP 14810 d herein marks 
the tenth system with three or more detected planets 
and only the sixth known to host three or more giant 
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(Msini > lOMggj^) planets. 4 The minimum masses 
and orbital periods of the planets in this systems are sim- 
ilar to those of v Andromcdae, but with the inner and 
outermost components reversed. For the three planets 
of HIP 14810 we find: M sini = 3.9, 1.2, and 0.6 M Jup , 

and P = 6.67, 147.7, and 952 d, respectively. We find 
modest but significant eccentricities for all 3 components 
(0.14, 0.16, and 0.17, respectively). 

Table [1] contains a summary of the stellar proper- 
ties of HIP 1481 (= BD+20 518), which sits at 53 pc 
(tt = 18.7 ± 1.3, Ivan Leeuwenl l2009h and has V = 8.5. 
We have performed an LTE analysis of our template 
spectra for HIP 14810 and derived its mass and ra- 
dius u sing the methods described in IValenti &: Fischer] 
(2005). Although HIP 14810 is a solar mass star (M = 
0.99M Q ), its metallicity ([Fe/H] = +0.26) and ev olution- 
ary st atus (AAfy = 0.63 mag, as calculated in [Wright 
(2004)) give it a spectral type of G5. Its low rotation 
(i>sini ~ 0.5 ± 0.5 km s _1 ) and Ca II H &K activity 
le vels (S = 0.16 , meas ured with the methods described 
in IWright et alJ (|2004T are consistent with it being an 
old star (age ~ 8 Gyr). This combined with its relatively 
small distance from the main sequence make it a partic- 
ularly good radial velocity targ et, since it is e xpected to 
exhibit very low levels of jitter (Wright 2005). 5 

3. VELOCITIES AND ORBITAL SOLUTION 

We began obser vations of HIP 1481 in 2005 as part of 
the N2K survey dFischer et al.l [2 005) at Keck Observa- 
tory using HIRES ( Vo et et apl994f ) and our usual iodine 
technique (|Butler et al.lll996T ) to achieve typical internal 
(random) errors of 0.8-1.4 ins -1 . The presence of the 
innermost planet and the large resuiduals to its orbital 
fit inspired the California Planet Search consortium to 
c ontinue regula r obser vations of this system at Keck. 

Wri ght et al.l (|2007l ) announced 6 the inner two plan- 
ets of HIP 14810, though with a rather poor fit for the 
c component due to the poor phase coverage and the 
unaccounted-for effects of the d component. Also compli- 
cating the fit was what is now obviously a spurious data 
point, acquired during early dusk when significant con- 
tamination from the Solar spectrum likely produced an 
erroneous radial velocity measurement. We have applied 
a more rigorous data retention scheme (based solely on 
the measured internal errors 7 , not deviations from a fit) 
to the data set presented in Table [2j These velocities and 
uncertainties supersede our previously published values 
for this star, as we continue to refine our data reduc- 
tion pipeline (| Wright et al .1 [2009h . Note that the times 

4 The others a re v Andomedae , 55 Cnc, u Ara e, HD 69830 
IILovis et al.ll2006T) . and HD 37124 llVogt et al.ll2005l) . 

5 We expect, based on hundreds of similar stars, 2 m s — 1 of 
jitter, consistent with the residuals to our best fit. This value i s 
somewhat lower than predicted by the formulae in IWright] (|2005T l 
because that work included instrumental sources of noise associated 
with the HIRES CCD detector in place prior to Aug 2004. The 
new detector has significantly better charge transfer properties, 
and apparently contributes a negligible amount to our overall error 
budget. 

6 A preliminary orbit for th e b component als o appears in the 
Catalog of Nearby Exoplanets (Butler et al. 2003). 

7 We derive radial velocity measurements from each of many 
independent "chunks" of spectrum for each observation. Our in- 
ternal error for a g iven observati on is determined from the variance 
of these velocities (Wright 2005). 



given in the table are in heliocentric Julian days, and the 
quoted errors are our internal (random) errors, with no 
"jitter" included. 

By 2007, residuals to a 2-planet fit clearly showed co- 
herent structure indicative of an outer companion. As 
Figure Q] shows, by late 2008, these residuals appeared 
to describe one complete orbit of a P ~ 950 d planet 
with modest eccentricity (e ~ 0.2). We have performed a 
Monte Carlo false alarm probability (FAP) analysis of the 
complete set of residuals to determine the likelihood that 
an orbital fit of this quality could have been arrived at by 
ch ance. This method i s very similar to the FAP analyses 
inlButler et ail (120061) fright et all (|2007fk iButler et all 
(|2009h : lHoward et all (|2009h . and we refer the reader to 
those works for more details. After binning the data in 
24- hour intervals and subtracting the best 2-planet fit, we 
redrew these residuals 1000 times (that is, we kept the 
times of observation the same for each trial, but at each 
time assigned a new velocity and velocity uncertainty 
pair randomly drawn from the entire set, with replace- 
ment). We then added the nominal orbital solution for 
the c component back into these new residuals and per- 
formed a thorough search for the best-fit 2-planet orbital 
solution to each of these 1000 realizations of the data. 
We measure the FAP as the fraction of these realizations 
for which we find a solution superior or equivalent in fit 
quality to the nominal solution. 

Note that this is an extremely conservative test of the 
FAP of the new planet, d, because we have not restricted 
the parameter search for the FAP trials to long period 
or low eccentricity planets. The best-fit solutions found 
for the artificial data sets are thus often at short peri- 
ods (P < 3d ) and/or high ecce ntricites e > 0.7 (see the 
discussion m IButler et all l2009). Nonetheless, even with 
this large parameter space available, we find that none 
of our 1000 realizations produced a better fit than the 
actual data, yielding an FAP < 0.1%. 

We have fit the data using the publicly available multi- 
planet RV-fitting I DL p ackage RV_FIT_MP, described in 
IWright fc Howard! (|2009l ). In Tabled we present the best 
3-planet Keplerian (kinematic) fit, which yields r.m.s. 
residuals of 2.3 ms _1 , and we plot the fit and velocities in 
Figure [T] We have observed just over one complete orbit 
for the outer component, HIP 14810d, and so its orbital 
parameters are sensitive to the assumption that there are 
no additional, external planets detectably influencing the 
velocities. In particular, it is possible that such an addi- 
tional planet could be contributing a nearly linear trend 
to the observed radial velocities, a trend which might 
be absorbed into the d component's orbital parameters 
as an inflated eccentricity. Any such degeneracy will be 
broken in the near future as the d component completes 
its second orbit. Continued observation of this system 
will thus reveal the presence of any additional detectable 
plan ets and allow for further analys is of their interactions 
(e.g. iFordl 12001 IWright et al.1l200l . 

4. DYNAMICAL MODELING AND STABILITY 
ANALYSIS 

We have performed long-term numerical integrations 
to test for stability of the orbital solution here, under 
the assumption that there are only three planets in the 
system. For these long-term stability tests, we applied 
direct n-body integrations to 100 orbital solutions consis- 
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Fig. 1. — Radial velocity curves for the HIP 14810 triple system. 




Planetary Orbital Phase 

Fig. 2. — Top Panel: The 1099 photometric observations of 
HIP 14810 in the combined Stromgren (b + y)/2 passband, ac- 
quired with the Til 0.8m APT over four observing seasons and 
plotted modulo the 6.673855-day orbital period of the inner planet 
HIP 14810b. Phase 0.0 corresponds to the predicted time of mid 
transit. A least-squares sine fit at the orbital period yields a semi- 
amplitude of only 0.00008 ± 0.00006 mag. Bottom Panel: The 
photometric observations of HD 14810 near the predicted time of 
transit replotted with an expanded scale on the abscissa. The solid 
curve shows the depth (0.01 mag) and duration (±0.0125 phase 
units) of a central transit, computed from the orbital elements 
and the planetary and stellar properties. The uncertainty in the 
predicted transit time is smaller that the diameter of the plotted 
points in the lower panel. The brightness levels of HIP 14810 within 
and outside the transit window agree to 0.00013 mag. Therefore, 
transits of HIP 14810b are ruled out for all reasonable densities. 



TABLE 1 
Stellar. Properties of HIP 14810 



Parameter 


Value 


Spectral Type 


G5 V 


RA 


hqIi 1 1 mi /iSoQn 

UO 14 4^±. ZOU 


Dec. 


+21 05 50. 49 


r>- V 


U. ( 


y 


8.52 


Distance (pc) 


52.9 ± 4.1 


M V 


4.9 


^cff(K) 


5485 ± 44 


logg [cms ] 


4.220 ± 0.06 


[Fe/H] 


+0.26 ± 0.03 


V sin i 


0.54 ± 0.5 km s" 1 


Mass (Mq) 


0.99 ± 0.04 


Radius (Rq) 


1.0 ± 0.06 


S 


0.16 


logRjjK 


-5.01 


AA/ V (mag) 


0.64 



tent with the RV data, assuming an edge-on orientation. 8 
We integrated for at least 10 s years u sing the hybrid in- 
tegrator in Mercury (| Chambers! 1 1 999f ) . For the majority 
of each integratio n, Mercury uses a mixed-v ariable sym- 
plectic integrator ( [Wisdom &: Holmanlll99li ) with a time 
step approximately equal to a hundredth of the Kcple- 
rian orbital period calculated at a semi-major axis equal 
to the pericenter distance of the closest planet. During 
close encounters, Mercury uses a Buhich-Stoer integra- 
tor with an accuracy parameter of 10 -10 . We identified 
each set of initial conditions as an unstable system if: 1) 
two planets collide, 2) a planet is accreted onto the star 
(astrocentric distance less than 0.005AU), or 3) a planet 
is ejected from the system (astrocentric distance exceeds 
100AU). 

All of our simulations proved stable, and we manually 
verified that in all cases the final orbits were qualita- 
tively similar to the initial conditions. We also ran an 
additional set of simulations for inclined (but coplanar) 
orientations of the system (i.e. smaller values of sini, 
and thus higher true planet masses and larger planet- 
planet interactions) and find the system to be stable for 
all tested scenarios with i > 3°. 

In addition, we performed dynamical analyses of three 
planet solutions using ensembles of initial conditions gen- 
erated via Bayesian posterior sampling methods (e.g. 
iFordl l2006f ) . both ignoring and including the planet- 
planet interactions. We found that the superposition 
of Keplerians approximation is a very good approxima- 
tion for the nominal edge-on configuration. As the orbits 
approach face-on, the planet-planet interactions eventu- 
ally become significant for inclinations of a few degrees. 
While we only considered coplanar configurations, we ex- 
pect that our results are likely representative for relative 
inclinations of up to 20 or 30 degrees. 



5. PHOTOMETRIC OBSERVATIONS 

We acquired 1099 good photometric observations of 
HIP 14810 during four observing seasons spanning 1180 

8 We have estimated the uncertainties in the orbital parame- 
ters in Table [3] using a variant of the bootstrapping described in 
I Butler ct al. (2006J). We generate a suite of plausible data sets and 
calculate the best-fit solution for each realization. We selected 100 
of these real izations as the basis for this stability analysis, as in 
IWrieht et aEI J2009T) . 
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TABLE 2 
Radial Velocities for HIP 14810 



HJD 


Velocity 


Uncertainty 


JD-2440000 


_ i 

m s 


_i 

m s 


1 3693 76627 


-147.5 


1.1 


1 3fi94 8371 


-489.3 


1 2 


1 369^ Q1 489 


-242.3 


I I 


1 3723 7908^ 


146.4 


1.0 


1 3794 (S9393 

J-O 1 ^'-t.UcfOUO 


304.8 


1 2 


1 374fi 81 748 


-13.1 


1.1 


1 3747 8^73 

J.O I *± I .OdJ I O 


-450.0 


9 


1 3748 73fi90 


-448.5 


1 2 


1 3749 741 8fi 


-86.0 


1.1 


1 37^1 990R8 
ioi oi. yuuuo 


342.3 


1 


1 3752 80978 


225.5 


0.8 


1 3753 f!9384 


-96.2 


I I 


1 3753 81 9fi1 

J.O I OO.OJ-^U J- 


-154.0 


1 


1 3753 9(1329 


-195.5 


1.2 


1 3775 83fi9fi 


-256.0 


9 


1 377(S 81 28(S 

lO 1 (U.O-LZiO'J 


107.0 


1.2 


1 3777 79391 

J.O i l i . t ^OUl 


332.0 


1 2 


1 3778 790R1 

J.O i 1 O. I iUUl 


398.1 


I I 


1 3779 7441 2 


220.8 


1.2 


1 39fi1 1 3098 


-293.8 




1 39fi2 1 3341 


-426.9 


1.0 


1 3981 971 79 


-492.7 


1 j 


1 "3989 94950 


-217.1 


1 2 


1 3983 98363 


137.6 


1.0 


1 3984 9991 5 
ioyo4:.uyyiu 


169.9 


1 2 


1 3984 98803 


333.0 


1.1 


1 3985 1 9437 


343.6 


1 2 


1 4993 9R949 
iiu^o. yuy^y 


88.2 


1 3 


1 4947 981 88 


-396.5 


1.2 


1 41 99 79fi95 


-208.9 


9 


1 41 "39 74873 

l-^kXOU . 1 '-to i o 


125.2 


9 


141 "31 84189 

14il01. Orfc J-0.ii 


334.8 


1 2 


1 41 38 89fi93 


336.4 


1 2 


1 4307 1 '3097 


92.8 


1.0 


1 433f3 9771 1 

IIOOU.UI 1 t -L 


-436.4 


1 


1 4344 031 62 


34.5 


1.0 


1 4345 1 4559 


306.4 


9 


1 439f3 89488 


-108.0 


1 


1 4397 91 1 1 2 


254.0 


1.0 


1 4398 8858fi 


411.3 


I j 


1 4399 89993 
iioyy.oy^yo 


364.3 


1 4 


1 4498 0901 3 


-184.5 


1 


14428.87298 


-470.7 


1.0 


14492.75897 


356.1 


0.9 


14544.74216 


260.4 


1.1 


14674.09015 


219.2 


1.2 


14718.05894 


146.8 


1.1 


14723.08578 


-434.2 


1.0 


14725.92426 


348.5 


1.2 


14727.05509 


271.3 


1.0 


14727.97892 


-45.6 


1.1 


14777.99836 


67.5 


1.1 


14805.83497 


358.4 


1.1 



days between 2005 November and 2009 February with 
the Til 0.8m automated photometric telescope (APT) 
at Fairborn Observatory. The Til APT and its two- 
channel Stromgren b and y photometer are very similar 
to the T8 0.8m APT and precision photometer described 
in iHenrvi (fl99l . 

The measurements of HIP 14810 were made differ- 
entially with respect to the comparison star HD 18404 
(V = 5.80, B — V = 0.42, F5 IV). We combined the 
Stromgren b and y differential magnitudes into a sin- 
gle (b + y)/2 passband to im prove the precision of each 
measurement. IHenrvi (fl999h gives further details on the 
acquisition, reduction, and calibration of the APT data. 

The standard deviation of a single observation from 



TABLE 3 

Orbital Elements for Exoplanets in the HIP 14810 System 



Property bed 



Per (d) 6.673855(19) 147.730(65) 952(15) 

T (JD-2440000) 13694.5980(70) 14672.2400(73) 14317.1980(73) 

e 0.14270(94) 0.164(12) 0.173(37) 

w (°) 159.32(38) 329.0(2.5) 286(19) 

K (ms" 1 ) 424.48(44) 50.01(46) 12.03(49) 

Msini (Mj ) 3.88(32) 1.28(10) 0.570(52) 

a (AU) 0.0692(40) 0.545(31) 1.89(11) 

r.m.s. (ms -1 ) 273" 

xl _ 1-01 
jitter (ms 1 ) 2 
^ohs §3 



Note. — For succinctness, we express uncertainties using par- 
enthetical notation, where the least significant digit of the uncer- 
tainty, in parentheses, and that of the quantity are to be under- 
stood to have the same place value. Thus, "0.100(20)" indicates 
"0.100 ±0.020", "1.0(2.0)" indicates "1.0 ±2.0", and "1(20)" indi- 
cates "1 ± 20" . 

the mean of the entire dataset is 0.00158 mag, which 
closely matches the typical measurement precision with 
this APT. Periodogram analysis of the full dataset finds 
no significant periodicities between one and several hun- 
dred days. In particular, least-squares sine fits to the 
6.673855- and 147.73-day orbital periods of the inner 
(b) and middle (c) planets yield semi-amplitudes of only 
0.00008 ± 0.00006 and 0.00021 ± 0.00006 mag, respec- 
tively, confirming that stellar activity is not the cause of 
the radial-velocity variations at these two periods. Al- 
though we do not expect radial velocity variations due to 
long period variatio ns in stellar activity in such old stars 
(Wrigl rFet~^[2008l) . we have also searched for photo- 
metric variations at the period of the outer (d) planet 
(P = 952 d). While our dataset is not yet long enough 
to make this determination definitively, we see no sug- 
gestion of photometric variations at this period so far. 

In the top panel of Figure [21 we plot the 
APT brightness measurements of HIP 14810 against 
phases computed from the 6.673855-day orbital pe- 
riod of the inner planet and a time of mid transit, 
JD 2,453,693.5856 ± 0.0022, predicted from the orbital 
elements in Table [3] The observations near the predicted 
time of transit are replotted with an expanded horizon- 
tal scale in the bottom panel of Figure [U The solid 
curve shows the depth (0.01 mag) and duration (±0.0125 
phase units) of a central transit of planet b, computed 
from the orbital elements and the planetary and stellar 
properties. The precisely determined orbit of HIP 148106 
(Figure[l] top panel) translates into an uncertainty in the 
mid-transit time that is smaller than the plotted points 
in the bottom of Figure^ The mean of the 120 observa- 
tions within the transit window is 2.81423±0.00015 mag; 
the mean of the 979 observations outside the window is 
2.81410 ± 0.00005 mag. Thus, the brightness levels in- 
side and outside of the predicted transit window agree to 
within 0.00013 mag, allowing us to rule out transits of 
the b component for all reasonable densities. 

6. CONCLUSIONS 

HIP 14810 is orbited by at least 3 giant planets, the 
outermost of which (HIP 14810 d) has completed just 
over one orbit since we began monitoring it at Keck Ob- 
servatory as part of the N2K project in 2005. This makes 



A Third Giant Planet Orbiting HIP 14810 



5 



it only the sixth system known to host more than 2 gi- 
ant (Msini > lOMggj^) planets. In retrospect, the 
previously published orbital solution for the middle (c) 
component was hampered by a spurious data point, poor 
phase coverage and the unaccounted-for effects of the d 
component. Continued monitoring will reveal the exis- 
tence of any additional outer planets. 

We have performed a dynamical analysis for a suite of 
orbital solutions assuming no fourth planet, and find the 
system to be stable at nearly any inclination (i < 3°). 
We have performed a photometric survey and find that 
the b component does not transit. 

The work herein is based on observations obtained at 
the W. M. Keck Observatory, which is operated jointly 
by the University of California and the California Insti- 
tute of Technology. The Keck Observatory was made 
possible by the generous financial support of the W.M. 
Keck Foundation. We wish to recognize and acknowl- 
edge the very significant cultural role and reverence that 
the summit of Mauna Kea has always had within the in- 
digenous Hawaiian community. We are most fortunate to 
have the opportunity to conduct observations from this 
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